Astronomy & Astrophysics manuscript no. 1.2mmG345.5 


©ESQ 2011 


November 1, 2011 





Massive star formation in the GMC G345.5+1.0: Spatial distribution 

of the dust emission 

Cristian Lopez 1 , Leonardo Bronfman 1 , Lars-Ake Nyman 2 , Jorge May 1 , and Guido Garay 1 

1 Departamento de Astronorma, Universidad de Chile, Casilla 36-D, Santiago, Chile, e-mail: clopez@das . uchile . cl 

2 European Southern Observatory, Casilla 19001, Santiago, Chile. 



o 

> 

O 

O 

ON 

(N 
< 

o 

Oh- 

6 

CO 



> 

m 
O 

in 



X 



ABSTRACT 

Context. Massive condensations in giant molecular clouds (GMCs) are linked to the formation of high mass stars, which are the 
principal source of heavy elements and UV radiation, playing an important role in the evolution of galaxies. 

Aims. We attempt to make a complete census of massive-star formation within all of GMC G345.5+1.0. This cloud is located one 
degree above the galactic plane and at 1.8 kpc from the Sun, thus there is little superposition of dust along the line-of-sight, minimiz- 
ing confusion effects in identifying individual clumps. 

Methods. We observed the 1.2 mm continuum emission across the whole GMC using the Swedish-ESO Submillimetre 
Telescope (SEST) Imaging Bolometer Array (SIMBA) mounted on the SEST. Observations have a spatial resolution of 0.2 pc and 
cover l°8x2°2 in the sky with a noise of 20mJy beam -1 . 

Results. We identify 201 clumps with diameters between 0.2 and 0.6pc, masses between 3.0 and 1.3xl0 3 M Q , and densities between 
5xl0 3 and 4xl0 5 cirr 3 . The total mass of the clumps is 1.2xl0 4 M Q , thus the efficiency in forming these clumps, estimated as the 
ratio of the total clump mass to the total GMC mass, is ~0.02. The clump mass distribution for masses between 10 and 10 3 M Q is 
well-fitted by a power law dN/dMocM - ", with a spectral mass index a of 1.7+0. 1 . Given their mass distribution, clumps do not appear 
to be the direct progenitors of single stars. Comparing the 1.2 mm continuum emission with infrared images taken by the Midcourse 
Space Experiment (MSX) and by the SPITZER satellite, we find that at least ~20% of the clumps are forming stars, and at most 
~80% are starless. Six massive-star forming regions (MSFRs) embedded in clumps and associated with IRAS point sources have 
mean densities of ~10 5 cm -3 , luminosities >10 3 L Q , and spectral energy distributions that can be modeled with two dust components 
at different mean temperatures of 28±5 and 200+10 K. 

Key words. ISM: clouds - ISM: stars: formation - ISM: dust, extinction - ISM: stars: circumstellar matter 



1. Introduction 

1.1. Giant molecular clouds 

High-mass stars are known to be born in massive and dense 
clumps embedded within giant molecular clouds (GMCs; 
Zinnecker & Yorke 2007). These GMCs have typical radii of 
~60pc, masses of 10 6 M and temperatures of ~10K, and are 
found to be concentrated toward the Galactic plane (Grabelsky 
et al. 1988). This gives rise to spatial and kinematic blending 
along the line-of-sight. For example, the CO emission, the most 
frequently used tracer of molecular gas for GMCs, with a low 
critical density (~ 10 2 cm -3 ), is often complex with multiple pro- 
file components. To overcome these observational difficulties 
and perform a complete census of massive-star forming regions 
(MSFRs) in a whole GMC, ideally requires one to study GMCs 
above and below the Galactic plane using a high density tracer. 
Recent surveys of dust condensations within whole GMCs have 
been made toward RCW 106 (Mookerjea et al. 2004), Cygnus 
X (Motte et al. 2007), and NGC 6334 (Munoz et al. 2007) using 
millimeter continuum emission as a high density tracer. 



1.2. Massive-star forming regions 

MSFRs are embedded in GMCs, generating a high Lyman con- 
tinuum photon flux (>2xl0 45 photons sec -1 ; Panagia 1973) that 
ionizes the surrounding gas and heats the surrounding dust. 
These MSFRs can thus be identified by infrared emission from 
heated dust and/or by radio emission from ionized gas. They can 



also be identified based on both their molecular lines from high 
density gas and their millimeter and sub-millimeter continuum 
emission from dust. 



Bronfman et al. (1996) showed that 60% of the IRAS 
point sources in the Galactic plane with far infrared (FIR) 
colors typical of ultra compact (UC) HII regions (Wood & 
Churchwell 1989) are associated with dense molecular struc- 
tures seen in the CS(2-1) line (>10 4 -10 5 cm -3 ). Hereafter, these 
sources will be referred to as IRAS-CS sources. From observa- 
tions in 1.2 mm continuum emission of 146 IRAS-CS sources, 
Faundez et al. (2004) showed that MSFRs are associated with 
condensations of gas and dust. Infrared studies, however, can- 
not provide a complete census of the birth sites of massive stars, 
since there are massive condensations that are undetected at in- 
frared wavelengths. For example, Garay et al. (2004) found four 
clumps with masses between 4xl0 2 and 2xlO 3 M and densi- 
ties of ~2xl0 5 cm -3 , without infrared emission, located close to 
clumps associated with MSFRs. They suggested that these cold 
(<17K), dense, and massive clumps will eventually form high 
mass stars. Hill et al. (2005) found 113 cold clumps, which have 
a mean mass of ~800M o , a mean radius of ~0.4pc, and a mean 
density of ~10 5 cm -3 . Beltran et al. (2006) found 95 cold clumps 
with a mean mass of 96 M Q , a mean radius of 0.4 pc, and a mean 
density of 9x1 4 cm -3 . 
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1.3. Clumps in GMCs 

Observations in molecular lines (e.g. Bains et al. 2006) and dust 
continuum emission (e.g. Munoz et al. 2007) on spatial scales 
smaller than ~ 1 pc show that GMCs have a fragmented struc- 
ture, and these sub-structures have been referred to as clumps 
(Williams et al. 2000). 

Clumps in GMCs have masses from 4 to 10 4 M , diameters 
from 0.2 to 2pc, and densities from 10 3 to 10 5 cirT 3 . The clump 
mass distribution is consistent with a power law dN/dMoc M - ", 
where dN/dM is the number of objects by mass interval, M is the 
mass, and a is the mass spectral index. The derived mass spectral 
indices range between 1.3 and 1.8 (e.g. Mookerjea et al. 2004; 
Munoz et al. 2007), values similar to those found for the mass 
distribution of molecular clouds as a whole (ar=1.5-1.6; Sanders 
et al. 1985; Solomon et al. 1987; Williams & McKee 1997). The 
similarity between the spectral mass indices suggests a common 
origin; however, different mechanisms have been proposed to 
explain the formation of clumps and GMCs. On the one hand, 
large-scale gravitational instabilities, in the combined medium 
of the collisionless stars and the collisional gas, drive spiral den- 
sity waves (e.g. Li et al. 2005), and are likely to be the main 
mechanism behind GMC formation. On the other hand, since 
GMCs are considered turbulent and supersonic, it is expected 
that the clump formation is, to first order, produced by ram pres- 
sures from supersonic flows, which can provide the seeds for a 
gravitational fragmentation (Ballesteros-Paredes et al. 2006 and 
2007; Bonnell et al. 2007; Klessen et al. 1998). 

1.4. This paper 

To undertake a complete census of dense and massive clumps, 
including those with and without infrared emission, we observed 
the whole of GMC G345. 5+1.0 in 1.2mm continuum emis- 
sion. The determination of physical properties of a sample of 
clumps belonging to a single GMC is desirable since, in this 
way, the distance to the GMC is a factor that influences neither 
the mass distribution or the relationships between physical prop- 
erties. However, these studies are difficult because of the super- 
position of dust and molecular structures along the line-of-sight 
to GMCs in the Galactic plane. We chose the GMC G345.5+1.0 
as our target for two reasons. First, it is located ~1° above the 
Galactic plane, so our observations in the 1.2 mm continuum 
emission and in the 13 CO(3-2) line are roughly free of confu- 
sion with background or foreground structures along the line- 
of-sight. Second, it is at 1 .8 kpc from the Sun, near enough to 
resolve clumps associated with MSFRs, but far enough away to 
permit a complete coverage of the GMC. 

The structure of this article is as follows: Sections [2] and [3] 
present, respectively, the main characteristics of G345.5+1.0 and 
the observations; Section 4 presents the results and a discussion, 
and Sect.[5]gives a summary of our study. 



2. GMC G345.5+1.0 

The GMC G345. 5+1.0 was first observed as part of the 
Columbia University - Universidad de Chile 12 CO(1-0) Survey 
of the Southern Galaxy (Bronfman et al. 1989). It is located 
approximately between 344?5 and 346?5 in Galactic longitude, 
and between 0°2 and 2°0 in Galactic latitude. Figure \T\ shows 
the spectrum of the 12 CO(1-0) line emission integrated over the 
area of the GMC. The emission of the GMC is between -33 and 
-2kms _1 (LSR velocities) with a peak at — 13.6kms _1 . 
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Fig. 1. Spectra of the 12 CO(1-0) line emission integrated over 
the whole area of the GMC G345.5+1.0 (between 344.°5 and 
346°5 in Galactic longitude and between 0.°2 and 2°0 in Galactic 
latitude; Bronfman et al. 1989). Emission from the GMC under 
study is between -33 and -2kms~' with a peak at -13.6 km s~'. 

Table 1. Summary of main characteristics of G345. 5+1.0. 



Distance 1.8 kpc 

Total mass 6.5x10 s M Q 

Radius 34 pc 

"Density 70 cm" 3 

"Column density 10 22 cm -2 

LSR velocity -13.6 km s -1 (between -33 and -2 km s -1 ) 



" The density and the column density are computed assuming a 
mean molecular weight of /j=2.29. 



We estimate the kinematic distance using the rotation curve 
determined by Alvarez et al. (1990), with a Galactocentric so- 
lar distance of 8.5 kpc and a solar LSR velocity of 220kms _1 . 
Considering a Galactic longitude of 345°5 and LSR velocity of 
-13.6kms _1 , the GMC is within the solar circle with two pos- 
sible kinematic distances: ~1.8 and 15 kpc. Thus, the GMC is 
~3 1 or 262 pc above the Galactic plane, a factor 0.5 or 4.4 of the 
HWHM of the molecular Galactic disk (~60pc; Bronfman et 
al. 2000), respectively. Therefore, 1.8 kpc is the most probable 
kinematic distance to the GMC. 

Physical properties of the GMC are estimated using the 
12 CO(1-0) line observations. The 12 CO(1-0) line emission in- 
tegrated over the full spatial and spectral extension of the 
GMC is 192 Kkms 1 deg 2 . Using a ratio of H2 column den- 
sity to integrated 12 CO(1-0) line emission NhJW co equal to 
1.56xl0 20 crrT 2 (Kkms 1 ) 1 (Hunter et al. 1997), the total mass 
of the GMC is 6.3xlO 5 M , corrected by a factor of 1.3 to ac- 
count for 25% of helium. Since its angular size is 2.°0xl°8, the 
GMC has a mean radius of ~34pc, mean column density of 
~10 22 cirT 2 and mean density of ~70crrT 3 , where the depth of 
the cloud has been assumed to be same as its radius. Table Q] 
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Galactic Long. 

Fig. 2. Map of the GMC G345.5+1.0 in 1.2 mm continuum emission. Observations were made using SIMBA, with a spatial resolu- 
tion of 0.2pc. They cover l°8x2°2 in the sky, with an rms of 20mJy beam -1 . 



summarises the main characteristics of this GMC. The derived 
physical properties confirm that the most probable distance to 
the GMC is ~1.8kpc. If at the far kinematic distance (15kpc), 
the GMC would have a mean radius of 280 pc and a total mass of 
4.4xl0 7 M o , values much higher than the typical value of 60 pc 
in radius and 1O 6 M in mass (Dame et al. 1986); Williams & 
McKee (1997) found that the molecular cloud mass distribution 
within the solar circle has an upper mass limit of 6xl0 6 M . 

3. Observations 

We observed the whole GMC G345. 5+1.0 in continuum emis- 
sion at 1.2 mm using SIMBA mounted on the SEST. The SEST 
is a 15-m diameter radio telescope, which operated between 70 



and 365 GHz, and SIMBA is a 37-channel hexagonal bolome- 
ter, operating at 250 GHz (1.2 mm), with a passband equivalent 
width of 90 GHz (FWHM). The configuration SEST-SIMBA had 
a beamsize of 24", which corresponds to a spatial resolution of 
-0.2 pc at the GMC distance (1.8 kpc). 

Observations were made in October 2002 and July 2003, us- 
ing the fast mapping mode, and consist of 185 images of 15' 
(azimuth) xlO' (elevation) in size. The scans were made in az- 
imuth at a rate of 80" s~', and they were separated in eleva- 
tion by 8"; the total integration time per map was about 25 
minutes. Measurements of the atmospheric opacity were made 
through sky dips about every three hours, and values at the zenith 
ranged between 0.09 and 0.31. Data were reduced using the 
MOPSI software (developed by Robert Zylka, IRAM, Grenoble, 
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Table 2. List of IRAS point sources along the line-of-sight to the 
GMC G345.5+1.0 observed in the CS(2-1) line by Bronfman et 
al. (1996). 



IRAS name 


Galactic coord. 


LSR velocity 1 


FWHM 1 




longitude 


latitude 


[kms -1 ] 


[kms -1 ] 


16533-4022 


344.845 


1.646 


undetected 


16571-4029 


345.208 


1.028 


-15.6 


4.9 


16577-4028 


345.286 


0.933 


undetected 


16575-4023 


345.332 


1.014 


-14.5 


2.8 


16561-4006 


345.393 


1.399 


-11.9 


3.5 


16557-4002 


345.395 


1.512 


-12.4 


3.1 


17009-4042 


345.490 


0.311 


-16.7 


6.0 


16562-3959 


345.494 


1.468 


-11.6 


5.5 


17008-4040 


345.499 


0.354 


-16.4 


4.8 


16596-4012 


345.717 


0.817 


-11.5 


4.9 



Values from Bronfman et al. (1996). 



France), and calibrated (in terms of flux density) with observa- 
tions toward Uranus, one made in October 2002, and an addi- 
tional eight made in July 2003. Bolometer channels were cor- 
rected for the correlated noise by inspecting the surrounding 
channels. The noise correlation between 1 and 900arcsec, and 
between 100 and 900 arcsec inferred an uncertainty/error smaller 
than 20% in the flux densities of detected sources. The calibra- 
tion factor has a value of 0.086 Jy counts -1 for October 2002, and 
a mean value of ~0. 069+0. 005 Jy counts -1 for July 2003. These 
variations agree with the uncertainty estimated by Faundez et al. 
(2004) of 20% in the flux density measurements using SIMBA 
data. 

For the 185 images, we achieved an rms of between 0.023 
and 0.080 Jy beam -1 , with a median of 0.036 Jy beam -1 , for a 
beam calibration area of ~653 arcsec 2 ("SIMBA Data Reduction 
Handbook"). The individual maps were combined in a mo- 
saic of L°8x2.°2 in size, centered at 345°40 in Galactic longi- 
tude and +1°10 in Galactic latitude, and with a final rms of 
~20mJybeam -1 . 



4. Results and discussion 

4. 1 . The 1 .2mm continuum emission 

Figure[2]presents the 1.2 mm emission image of the whole GMC 
G345.5 + 1.0. The total flux density of the GMC G345.5+1.0 is 
~365 Jy, which is estimated by integrating the intensity over the 
whole area of the cloud. 

As can be seen in the integrated spectrum of the 12 CO(l - 0) 
line (Fig. [TJ, across GMC G345. 5+1.0 (LSR velocity between 
-33 and -2kms -1 ), there are additional molecular gas compo- 
nents along the line-of-sight, particularly in the ranges -170 to 
-lOOkrns -1 , -75 to -33kms -1 , and to 5kms -1 . Hence, the 
question arises as to whether the 1.2 mm continuum emission 
only traces dust condensations within GMC G345.5+1.0? To ex- 
amine the association of the GMC with 1 .2 mm continuum emis- 
sion, Fig.[3]shows images of the velocity-integrated 12 CO(l - 0) 
emission in three velocity ranges (-200 to -33 kms -1 , -33 to 
-2 kms -1 , and -2 to 50 kms -1 ) superimposed with contours of 
the 1 .2 mm continuum emission. Figure[3]shows that most of the 
emission detected in 1.2 mm is associated with the GMC. About 
1% of the total observed area might also be associated with gas at 
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Fig.3. Integrated 12 CO(1-0) emission toward GMC G345.5+1.0 
in different LSR velocity ranges (Bronfman et al. 1989). Top: 
from -200 to -33 kms -1 . Middle: from -33 to -2 kms -1 . 
Bottom: from -2 to 50 kms -1 . Magenta circles mark spatial 
and spectral positions of detections in the CS(2-1) line toward 
MSFRs (Table 0. Contours represent 1.2 mm continuum emis- 
sion at 5 times rms, ~0.1 Jy beam -1 . 
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Fig. 4. IRAS point sources along the line-of-sight of the GMC 
G345.5 + 1.0 observed in the CS(2- 1 ) line (Bronfman et al. 1996). 
Gray scale represents 1 .2 mm continuum emission. Arrows mark 
CS(2-1) line observations, and crosses indicate observations 
without detection (see Table|2]). 



velocities <-33kms 1 (Fig.[3j top), localized mainly in the re- 
gion of the IRAS point sources 17008-4040 (G345.499+0.354) 
and 17009-4042 (G345.490+0.31 1). 

From the survey of Bronfman et al. (1996), we find that 
there are eight IRAS-CS sources within the region and two 
IRAS point sources that are not detected in the CS(2-1) line (see 
Table|2]i. As is shown in both the 12 CO(1-0) maps (Fig. [3} and 
the 1.2 mm continuum emission map (Fig. |4), these MSFRs, 
or IRAS-CS sources, are associated with the GMC and have a 
counterpart in 1 .2 mm. They correspond to the most dense and 
massive dust condensations (see Table 0. The two IRAS point 
sources not detected in the CS(2-1) line were also not detected 
in the continuum (see Fig. The eight IRAS-CS sources 
include the IRAS point sources 17008-4040 and 17009-4042. 
Line profiles toward these two objects in the 12 CO(1-0) and 
CS(2-1) lines are shown in Fig. [5] Gas components with 
velocities <-33kms _1 observed in the 12 CO(1-0) line are not 
observed in the CS(2-1) line, suggesting that they correspond to 
regions of low density gas. Since 1.2 mm continuum emission 
traces high densities (e.g. Faundez et al. 2004), it should not 
be detected in these clouds. In summary, from observations in 
the 12 CO(1-0) and CS(2-1) lines, we conclude that the 1.2 mm 
continuum emission is associated only with the GMC. 
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Fig. 5. Line profiles toward the IRAS point sources 17008-4040 
and 17009-4042: top image, 12 CO(1-0) line profiles (Bronfman 
et al. 1989), over a map of their integrated emissions (color 
scale) and with contours of 1 .2 mm continuum emission; middle 
and bottom images, CS(2-1) line profiles (their observing posi- 
tions are indicated as black dots in the top image; Bronfman et 
al. 1996). 



4.2. Identification of clumps 

The structure of the GMC observed in 1 .2 mm continuum emis- 
sion (Fig. [2]) is fragmented, and it is possible to distinguish sev- 
eral clumps. 



To identify clumps we utilize CLUMPFINEQ (Williams et 
al. 1994), which creates contours over data, searches for peaks of 
emission to locate clumps, and follows them down to the lower 
intensity contour. 

1 http://www.ifa.hawaii.edu/users/jpw/clumpfind.shtml 
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CLUMPFIND finds 201 clumps in the 1.2 mm continuum 
emission map of the GMC, containing ~100% of the total emis- 
sion above 3<x. We used a lower intensity contour of three rms, 
~0.06 JybeatrT 1 , and a contouring interval equal to twice the 
rms noise, ~0.04 JybeatrT 1 . To delete fictitious structures, we 
imposed two conditions on the CLUMPFIND output, that the 
angular size of the emission and the emission peak of clumps 
had to be greater than the beam size, ~24"x24", and five times 
rms, ~0.1 Jy beam" 1 , respectively. The angular size is defined 
to be the angular area inside the lowest intensity contour (three 
rms). The 201 identified clumps have areas between ~0.18 and 
7.3arcmin 2 , emission peaks between 0.1 and 9 JybeatrT 1 , and 
flux densities between 0.089 and 40 Jy. 

Table[5] available online in electronic form at the publisher^ 
(EDP Sciences), shows the characteristics of each clump calcu- 
lated in this section, Sect. 14.31 and Sect. 14.41 Column 1 gives 
clump names; columns 2 and 3, Galactic coordinates of peaks 
in 1.2 mm continuum emission; column 4, 1.2 mm flux densi- 
ties; column 5, diameters (deconvolved FWHM sizes); column 
6, masses; column 7, densities; column 8, column densities; and 
column 9, if clumps have ("Y") or do not have ("N") an infrared 
counterpart from MSX and SPITZER observations. 

4.3. Physical properties of clumps 

First, we estimate the minimum gas column density that can be 
detected given the rms noise of our observations. Assuming that 
1 .2 mm continuum emission is optically thin and produced by 
dust, the column density N is (Hildebrand 1983) 



1000 



N 



S 1.2 mm Rgd 



Clfi m H k\2mm Bl.2mm(Tdust) 



(1) 



where Q is the beam solid angle, Si. 2mm is the flux density at 
1.2 mm, fi is the mean mass per particle, equal to ~2.29 for an 
H2 cloud with a 25% contribution of helium (Evans 1999), m# 
is the hydrogen atom mass, ki.2 mm is the dust absorption co- 
efficient at 1.2mm, equal to ~lcm 2 g~' for protostellar cores 
(Ossenkopf & Henning 1994), B\2mm{Tdust) is the Planck func- 
tion at both 1 .2 mm and a dust temperature Td usl , equal to ~30 K 
for regions of massive-star formation (Faundez et al. 2004), and 
R g d is the ratio of gas to dust mass, ~100 (Hildebrand 1983). For 
a solid angle limit of 24"x24" and an intensity limit of five rms, 
~0.1 Jy beam" 1 , the minimum flux density is ~88mJy and the 
minimum column density that can be detected is ~4xl0 21 cm -2 , 
which corresponds to a visual extinction of 4 mag, assuming that 
N H2 /A v ~10 2l cm- 2 mag~ l (Bohlinetal. 1978). 
Masses of clumps, M c , are estimated as 



M.-j 



1.2 mm "grf 



Red d~ 



Bl.2mm(Tdust) 2 I. 



where dA is the differential element of area (dA= d 2 dQ) and 
d is the distance to the GMC (~1.8kpc). Since we insist that 
identified clumps have intensities and dimensions greater than 
0.1 Jy beam -1 and 24"x24" respectively, the lower limit to their 
masses is ~2.9M . The derived masses of clumps range from 
3.0 to 1.3xl0 3 M o , with a total mass of 1.2xl0 4 M Q (see Tables 
[3] and [5]). The efficiency in forming these clumps, estimated as 
the ratio of the total clump mass to the total GMC mass, is thus 
-0.02. 

The clump mass distribution (CMD) is shown in Fig.|6l plot- 
ted as dN/dlog(M/M e ) versus mass, where dN/dlog(M/M e ) 



Fitted a = 1.7±0.1 
IMF a-2.35 
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Fig. 6. Mass distribution of identified clumps in G345.5+1.0, 
plotted as dN/dlog(M/M ) versus mass, where dN/dlog(M/M ) 
is approximated by the number of clumps AN within a log- 
arithmic mass interval Alog(M/M ). Here, Alog(M/M ) is 
constant,~0.44. Error bars are estimated by sjAN/A\og(M/M Q ). 
The arrow shows the clump mass limit, ~2.9M . The 
continuous line represents the mass distribution fit with 
dN /dlog(M /M Q )ccM l - a , where the spectral mass index a is 
1.7+0.1 for masses between ~10 and 1.3xlO 3 M . The dashed 
line displays the spectral mass index for the stellar initial mass 
function (IMF) of the solar neighborhood for stellar masses 
greater than 0.5 M (e.g. Kroupa 2002); the line is forced to pass 
through the peak of the clump mass distribution. 



is approximated by the number of clumps, AN, within a loga- 
rithmic mass interval Alog(M/M ). In this figure, Alog(M/M ) 
is constant, at a value of ~0.44. The CMD is well-fitted by a 
power law dN/dlog(M) ocM _a+1 , which can be expressed as 
dN/dM ocM~'\ with the spectral mass index, a, equal to 1.7+0.1 
for masses between ~10 and 1.3xlO 3 M . The correlation coef- 
ficient of the fit is 0.993. Between 10 M Q and 1.3xlO 3 M , bin 
size variations of Alog(M/M ) between 0.18 and 1.1 result in 
values of a consistent with 1.7+0.1. Since a is 1.7, the popula- 
tion is dominated by clumps with low masses, but the total mass 
is dominated by the most massive clumps; for example, 50% of 
the population is between 10 and 27 M , but contains only 10% 
of the total mass. The turnover below ~1OM is produced by 
an incompleteness of the clump catalog caused by the combina- 
tion of the spatial resolution and flux density limit of the survey. 
However, observations of higher spatial resolution (<0.01pc) 
could result in a spectral mass index of ~2.35, resolving core 
structures (e.g. Motte et al. 1998). Beltran et al. (2006) studied a 
sample of IRAS sources associated with MSFRs in 1 .2 mm con- 
tinuum and found a spectral mass index of 1 .5 for clumps with 
masses between ~1OM and 1O 2 M and 2.1 for clumps with 
masses between ~10 2 and 1O 4 M . For clumps identified here 
with masses higher than 100 M , we find a spectral mass index 
of 1.6+0.1, as is shown in Fig. [7] in agreement with the previous 
fit considering clumps with masses between 10 and 1.3xl0 3 M . 
It is necessary to observe more whole GMCs to confirm these re- 
sults. 

Clump diameters, D c , are estimated from the deconvolved 
FWHM size of their emissions. We used the FWHM size Qfwhu 
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Table 3. Summary of the physical properties of the identified 
clumps. 



Mass [M ] 



Fig. 7. Mass distribution of identified clumps in G345. 5+1.0 
with masses higher than 100 M , plotted as dN/dlog(M/M Q ) ver- 
sus mass, where dN/dlog(M/M ) is approximated by the number 
of clumps AN within a logarithmic mass interval Alog(M/M ). 
Here, Alog(M/M ) is constant,~0.37. Error bars are estimated 
by -J AN I A log(M/M ). The continuous line represents the mass 
distribution fit with dN/dlog(M/M Q )ocM l ~ a , where the spec- 
tral mass index a is 1.6+0.1 for masses between ~100 and 
1.3xlO 3 M . 

estimated by CLUMPFIND algorithm, thus 



where % eam is the beam-size. Considering a distance of 1.8 kpc 
to the GMC and clumps that have a reliable D c , i.e. D c >6b eam , 
clumps have diameters between 0.2 and 0.6 pc. 

From the masses and sizes, and assuming a spherical and ho- 
mogeneous density distribution, we estimate mean clump densi- 
ties, using the expression 

p — Id niH n, 

where p is the mass density and n is the particle density. 
Densities of clumps are between 5xl0 3 and 4xl0 5 cirT 3 . Mean 
column densities of clumps, N c , are estimated as 



N, 



M c 



lim H n(DJ2) 2, 



and range between 4xl0 21 and 4xl0 23 °™ 2 



' cm Tables [3] and [5] 
show physical properties for each clump and a summary of them, 
respectively. 

Figure[8]shows a plot of mass versus diameter for the clumps 
detected toward GMC G345.5+1.0 with reliable diameters. The 
dotted lines indicate constant densities at 10 3 , 10 4 , 10 s , and 
10 6 cirT 3 . The majority of clumps have densities between 10 4 
and 10 5 cm 4 . 

The physical properties of detected clumps are similar to 
those found in other GMCs (e.g. Mookerjea et al. 2004). 



4.4. Association with infrared emission 
(IRAS-MSX-SPITZER) 

Stars form in clumps, heating their surrounding dust, which re- 
radiates at infrared wavelengths. This is illustrated in Fig.[9]that 
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shows a strong spatial correlation between the 1.2 mm contin- 
uum emission and infrared emission at 21.34yum (from MSX 
observations). To quantify the correlation, we searched for in- 
frared emission inside clump emission areas, using MSXB im- 
ages at 8.28, 12.13, 14.65, and 21.34/mi and SPITZEB0 (IRAQ 
images at 3.6, 4.5, 5.8, and 8.0/an. We find that -20% of all 
clumps have an infrared counterpart in all MSX and SPITZER 
bands (see Table |5j. The rest of the clumps, ~80%, are not de- 
tected in all MSX and SPITZER bands, particularly not in 12. 13, 
14.65, and 21.34/an. Since 8.0/mi MSX band and SPITZER 
IRAC bands are sensitive to the polycyclic aromatic hydrocar- 
bon (PAH) emission and the photospheric emission from stars 
(e.g. Chavarrfa et al. 2008), clump not detected in all MSX and 
SPITZER bands are considered to have no counterpart at in- 
frared wavelengths. Since both the MSX and SPITZER bands 
have sensitivity limits, the percentage of detections is a lower 
limit to the number of clumps that are forming stars, and the 
percentage of failed detections is an upper limit to the number of 
clumps that are not forming stars. 

Nine clumps are associated with six IRAS point sources 
classified as MSFRs with luminosities >1O 3 L (see Sect. 14.51 
and Table @). As Fig. [8] shows, these clumps have densi- 
ties of ~10 5 cm -3 , suggesting that there is a threshold density 
above which massive stars can form. These values are consis- 
tent with the typical density of clumps associated with MSFRs 
(~10 5 cm 4 ; Faundez et al. 2004). 

As Fig.[8]shows, clumps that emit detectable infrared emis- 
sion tend to be more massive than remaining clumps. Clumps 
without infrared emission (cold or starless clumps) have a mean 
mass of 21 M , and clumps with an infrared counterpart have a 
mean mass of 2. 1 x 10 2 M . Furthermore, all clumps with masses 
higher than ~200 M have an infrared counterpart. 

MSX point sources associated with clumps within the GMC 
G345. 5+1.0 have mid-infrared colors S2\/S& from 0.9 to 30, 
Si4/Sx from 0.4 to 8, and S 12/Sg from 0.7 to 4, where Sg, 
S 12, S 14, and S21 are the flux densities at 8.28, 12.13, 14.65, 
and 21.34 /im, respectively. These ratios cover those of massive 
young stellar objects (MYSOs; Lumsden et al. 2002): S 21/^8 > 
2 and S 21 > S 14 > S%. About 7% of the clumps contain MSX 
sources that satisfy this criterion, and these clumps have masses 
>36M (see Fig. EJ. 

The existence of clumps with and without infrared emission 
suggests that clumps in the GMC are at different evolutionary 
stages. Cold clumps have masses between 3.0 and 1.9xlO 2 M , 
where the most massive ones are possible progenitors of MSFRs. 
For example, we estimate that the least massive clump associated 
with a MYSO has a mass of ~36 M , and we identify seven cold 
clumps with densities >10 5 cirT 3 and masses >36M , which 
will eventually collapse to form high-mass stars. 



http://irsa.ipac.caltech.edu/ 
http://irsa.ipac.caltech.edu/ 
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Fig. 8. Mass versus diameter for the clumps detected toward the GMC G345. 5+1.0 in 1.2mm continuum emission with reliable 
diameters. Filled circles indicate clumps detected in infrared MSX and SPITZER bands. Open circles indicate clumps that do not 
have an infrared counterpart. Triangles indicate clumps associated with MSFR-IRAS sources, which have luminosities >10 3 L o . 
Boxes indicate clumps associated with MSX sources that satisfy MYSO candidate criterion (Lumsden et al. 2002). Arrows mark 
detection limits for masses (~2.9M ) and diameters (-0.2 pc). The continuous line indicates the detectable mass as a function of 
diameter (sensitivity limit), considering an intensity limit of five rms (~0.1 Jy beam -1 ). Dotted lines indicate mean densities at 10 3 , 
10 4 , 10 5 and 10 6 cm 4 . The densities are computed assuming a mean molecular weight of yi-2.29. 



Do clumps form single stars? One way to assess this is to 
compare the slope of the clump mass distribution with that of 
the stellar initial mass function (IMF) (e.g. Motte et al. 1998; 
Lada et al. 2007). Equal slopes would indicate that the origin of 
the stellar IMF has its direct roots in the origin of the clump mass 
distribution. The spectral mass index a of the clump mass dis- 
tribution determined here is consistent with that of other investi- 
gations (e.g. Munoz et al. 2007), but differs from that estimated 
for the stellar IMF of the solar neighborhood for stellar masses 
higher than 0.5 M {a~235\ e.g. Kroupa 2002). This suggests 
that the detected clumps do not directly form stars, and other 
processes are necessary to determine the stellar initial masses, 
such as the fragmentation of clumps, mainly of the most mas- 
sive ones. Figure [6] compares the IMF spectral mass index with 
the clump mass distribution. 



4.5. Dust properties of massive-star forming regions 
associated with clumps and IRAS point sources 

Regions of massive star formation are embedded in massive 
clumps, and the intense Lyman flux produced by them heats the 
surrounding dust, which re-emits mainly at far infrared wave- 
lengths with characteristic colors. To study the physical proper- 
ties of dust in these regions, we examine the spectral energy dis- 
tributions (SEDs) of MSFRs associated with IRAS point sources 
and clumps detected here, assuming that their emissions are from 
dust. 

Within the GMC, there are eight MSFRs associated with 
IRAS-CS sources (Table 0. We added one more source, IRAS 
16533-4009, which is embedded in 1.2 mm continuum emis- 
sion, has a high luminosity, ~9xl0 4 L , and increasing IRAS 
flux densities from 12 to 100/jm; however it does not satisfy 
the far-infrared color criterion defined by Wood & Churchwell 
(1989), since its flux density in 25 fim is an upper limit. Figure 
[TT] shows 8.0/iin images from SPITZER data with contours 
of 1.2 mm continuum emission for all these sources. For six 
of these objects, the SPITZER emission is embedded within 
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Fig. 9. Image in 2 1 .34 jjm from MSX observations toward GMC 
G345.5+1.0 with contours of 1.2 mm continuum emission at 
three times rms, ~0.06 Jy beam -1 . 



1 .2 mm continuum emission. For more reliable estimates in our 
SED study, we only consider these six sources: IRAS 16533- 
4009, IRAS 16562-3959, IRAS 16571-4029, IRAS 16596-4012, 
IRAS 17008-4040, and IRAS 17009-4042. 

Figure[T2]displays the six SEDs constructed using our obser- 
vations in 1.2 mm continuum emission, infrared data at 12, 25, 
60, and lOOjurn from the IRAS Point Source Catalog (version 
2.0), and at 8.3, 12.1, 14.7 and 21.3 /an from the MSX Point 
Source Catalog (version 2.3). Because observations were per- 
formed using different beam sizes, for IRAS ~300", for MSX 
~20", and for 1.2 mm 24", we consider all emission within the 
IRAS beam; thus, two SEDs are associated with more than one 
1.2 mm clump. 

The SEDs of MSFRs can be modeled as several dust compo- 
nents at different temperatures (e.g. Faundez et al. 2004; Morales 
et al. 2009). Because of the shape of the six SEDs, we model 
them as two dust components at different temperatures, cold and 
warm components, including the absorption of the radiation by 
assuming that the warm component is embedded in the cold one. 
The total flux density, S at frequency v is approximated by 

c total c cold , c warm 

J v ~ O v T & v ) 



where 



s coM = a cold By(T cold )(l _ eX p(- T » W )) 



and 



10000 
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- 
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Fig. 10. Dust opacity spectrum utilized in the SED models. It 
was estimated by Ossenkopf & Henning (1994)", assuming a 
Mathis-Rumpl-Nordsieck initial size distribution with thin ice 
mantles and 10 5 yr of coagulation at a gas density of 10 5 cm 4 . 
a |http://vizier,u-s trasbg.fr/viz-bin/VizieR| J/A+A/291/943/tablel 



The parameters S c v om , Q COM , T c d °f ( , and -f° M are the flux den- 
sity, the solid angle, the dust temperature, and the optical depth 
of the cold component, respectively, and Sy arm , €l warm , T™™\ 
and t™ 1 '" 1 are the flux density, the solid angle, the dust tempera- 
ture, and the optical depth of the warm component, respectively. 
In addition, B v {T c d °^ t ) and B Y (T™™) are the Planck function at 
dust temperatures T c . old and T d '"™, respectively. For both com- 
ponents, £1 can be expressed as 

Q = 7r (9/2) 2 , 

where 8 is the angular diameter. The optical depths are given by 
(e.g. Evans 1999) 



cold 



N coU , 



and 



N. 



dust 



where N c f w , and m a ™ 

dust dust 



are the dust column densities in g cm -2 
for the cold and warm components, and k y is the dust opac- 
ity. We use dust opacities estimated by Ossenkopf & Henning 
(1994) for protostellar cores. They computed opacities consid- 
ering the Mathis-Rumpl-Nordsieck (MRN) distribution for the 
diffuse interstellar medium (Draine & Lee 1984) as the initial 
size distribution for dust, without and with ice (thin and thick), 
and without and with coagulation (after 10 s years for densities 
between 10 5 -10 8 cm -3 ). In the case of regions with ice depletion 
produced by the heating of central sources, they recommended 
opacities for the model with thin ice mantles and coagulation for 
a density of 10 5 cm -3 . These dust opacities are shown in Fig.fTOl 
for frequencies between ~2.3xlO n and 10 14 Hz. 

Thus, in our SED model, each dust component has three val- 
ues to fit of Td USt , 9, and Ndust- However, the dust column density 
of the warm component, N d v "' s "\ is difficult to estimate, because 
it is more sensitive to the emission in the Rayleigh- Jeans part of 
the spectrum (hv«kT warm ), where the emission is dominated 
by the cold component. To overcome this problem, we assume 
that the two components have equal densities, thus 



N. 



dust 



/warm 
U di(st 
rcold 
V dusi 



cold 



c warm c\warm n /< r rwarm\/A ^warm\\ ^.cold /<->\ 

S y =12 By(T dust )(1 - exp(~T„ ))exp(-r v /2). 



Given the simplicity of the SED model and the poor sensitiv- 
ity of the data to A^°™, a more realistic density distribution is 
unnecessary. 
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To enable a more reliable comparison, angular diameters are 
converted into spatial diameters, and dust column densities to 
gas column densities. Thus 



T\cold _ J (fO 
u dust - " v dust' 



cold 



D 



fiat 
dust 



and 



n: 



cold 



j swarm 
d Vdust ' 



jim H 



where (f, ok [ and Q wa ™ are in radians. In this way, our model has 

dust dust J ' 

five variables: dust temperatures and sizes for the two compo- 
nents, and gas column density for the cold component. 

Tableland Fig.[T2]display the results of the fits. The mean 
dust temperatures of each component are 28+5 K (cold) and 
200+ 10 K (warm). The sizes and column densities of the cold 
component agree with those estimated by 1.2 mm continuum: 
sizes vary by a factor of 0.7-1.5 and column densities vary by a 
factor of 0.5-3. Estimates of luminosities, from the integration of 
SED models, are >10 3 L . Given the sizes and column densities, 
the total mass is dominated by the cold component (~99% of 
the total mass), and is similar to that estimated from the 1 .2 mm 
continuum emission, varying by a factor of 0.8-1.6. 

Dust characteristics of clumps associated with MSFRs es- 
timated in this paper are consistent with previous works (e.g. 
Faundez et al. 2004; Molinari et al. 2000; Molinari et al. 2008), 
where the cold dust temperature in regions of massive star for- 
mation is ~30 K. 

The SED models have a discrepancy with data in ~12//m 
bands (see Fig. [121) . which can be produced by not consider- 
ing the PAH emission in the SED models (e.g. van Dishoeck 
2004), or by an excess in the dust opacities utilized at these 
wavelengths, affecting the modeled radiation from the hot com- 
ponent that is absorbed by the cold one. 

The dependence of SED models on variations in the fitted pa- 
rameters is shown in Fig.[T3j which displays the SED for clump 
1, or IRAS 16562-3959, with the best-fit model and models with 
variations in the fitted parameters Tf u „ D c f ld t , 7T 1 ™ W. ar f, and 

1 dust' dust' dust ' dust 

N c °Jf. For each dust component, variations in its temperature af- 
fect both its peak emission frequency and its luminosity, both 
values increasing as the temperature increases. Variations in ei- 
ther or N™'f affect the luminosities of the two components; 
when D c ^jj t or N c °Jf increases, the luminosity of the cold dust 
component increases, whereas radiation absorption of the warm 
component also increases, reducing its luminosity. When D^"™ 
increases, the luminosity also increases. 



5. Conclusions 

We have robustly detected the whole of the GMC G345.5+1.0 in 
1 .2 mm continuum emission at a spatial resolution of 0.2 pc, and 
conclude that: 

- The GMC is fragmented. We have identified 201 clumps, 
which have beam-corrected diameters between 0.2 and 
0.6 pc, masses between 3.0 and 1.3xlO 3 M , and densities 
between 5 x 1 3 and 4x 1 5 cuT 3 . 

- The total mass of the clumps is ~1.2xl0 4 M Q , and after 
comparing with the total mass of the GMC of ~6.5xl0 5 M , 
we inferred that the efficiency in forming these clumps is 
-0.02. 



The clump mass distribution is well-fitted by a power law 
dN/dMocM~ ff , where the spectral mass index a is 1.7+0.1. 
The total mass is dominated by massive clumps, but the 
population is dominated by clumps with low masses. 

The spectral mass index of the clump mass distribution is 
different from that of the stellar IMF. Thus our detected 
clumps are probably not the direct progenitors of single stars. 

Comparing with MSX and SPITZER (IRAC -bands) obser- 
vations, 20% of the clumps have an infrared counterpart 
in all MSX and SPITZER bands. The remaining clumps, 
~80%, are considered to have no counterpart at infrared 
wavelengths. The percentage of detection is a lower limit to 
the number of clumps forming stars, while the percentage of 
no detections is an upper limit to the number of clumps that 
are not forming stars. 

Regions of massive-star formation within the cloud, as- 
sociated with IRAS point sources, have SEDs that can 
be modeled with two dust components at different mean 
temperatures of 28±5 and 200+ 10 K. 
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Fig. 11. Images of 8.0jum emission (SPITZER data) toward clumps detected in 1.2 mm continuum emission and associated with 
IRAS point sources. Contours represent 1.2 mm continuum emission at 0.06, 0.12, 0.24, and 0.48 Jy beam -1 (rms is 0.02 Jy beam -1 ). 
IRAS source names are given at the top of each image, and clump numbers are indicated at the peak of 1 .2 mm continuum emission. 
Red circles are centered on the coordinates of IRAS point sources, with diameters of 5' (an approximation of the angular resolution 
of IRAS observations at lOOjum). 
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Fig. 12. The SEDs of massive-star forming regions associated with massive clumps detected in 1.2 mm continuum emission; top 
labels show names of the clumps. Dots with error bars are flux densities estimated from SIMBA, IRAS, and MSX observations. 
Each SED is modeled with two dust components at different temperatures (physical parameters for each model are in Table |4); 
drawn lines are the total flux density of the two dust components, and dashed lines are the contributions of each dust component. 
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Table 4. SED models for MSFRs associated with IRAS point sources and mas- 
sive clumps detected in 1.2 mm continuum emission. Each model consists of two 
dust components with equal densities at different temperatures (cold and warm 
components). Column 1 shows names of IRAS point sources; Column 2, names 
of clumps; in Columns 3 to 8, we show the fitted physical parameters: dust tem- 
perature, diameter, and column density, respectively, for the cold (Cols. 3-5) and 
warm (Cols. 6-8) components; Column 9, masses of the cold component; and 
Column 10, total luminosities. 



IRAS 


Clumps 


Tdust 

K 


Cold 
pc 


N 

^ yets 

10 23 crrT 2 


Tdust 

K 


Warm 
pc 


N 

ly gas 

10 22 cm" 2 


Cold mass 
10 3 M o 


Luminosity 
1O 4 L 


16562-3959 


1 


32 


0.67 


1.6 


194 


0.03 


0.69 


1.0 


5.3 


17009-4042 


2 


29 


0.86 


1.6 


202 


0.03 


0.51 


1.7 


4.7 


16533-4009 


3,7, 10 


25 


0.79 


3.2 


205 


0.07 


2.7 


2.8 


9.4 


17008-4040 


4 


28 


0.68 


1.5 


205 


0.03 


0.57 


1.0 


2.9 


16571-4029 


5,9 


34 


0.58 


1.3 


209 


0.03 


0.65 


0.64 


5.3 


16596-4012 


12 


22 


0.34 


2.2 


203 


0.008 


0.49 


0.36 


0.22 
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Variations in the cold temperature 



Variations in the cold diameter 
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Fig. 13. Dependence of the SED model on variations in the fitted parameters. Plots show the SED for clump 1 with the best-fit 
model (see Table |4]i and variations in each parameter, for T%££, D c J^ t , T^"™', D™™' and Ng°'f, when increasing and decreasing the 
best-fit value. 
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Table 5. Properties of identified clumps in the GMC G345.5+1.0. Column 1 
gives clump names; columns 2 and 3, Galactic coordinates of peaks in 1.2 mm 
continuum emission; column 4, 1.2 mm flux densities; column 5, diameters; col- 
umn 6, masses; column 7, densities; column 8, column densities; and column 9, 
if clumps are detected in all infrared MSX and SPITZER-IRAC bands ("Y") or 
not ("N"). The densities and the column densities are computed assuming a mean 
molecular weight of /j=2.29. 



Name 


Galactic coord. 


S 1 .2mm 


Diameter 


Mass 


Density 


Column density 

_9 

cm - 


Infrared 




long. 


lat. 


Jy 


pc 


M Q 


cm 


counterpart 


1 


345.490 


1.471 


\ a 

40 


a c 

0.5 


1 1 •. * 1 a3 

1.3x10 


'"l A\ s 1 A*i 

2.9x10 


i a23 
3.2x10 


Y 


2 


345.488 


0.316 


1 o 

.18 


a c 

0.6 


1 1 - 1 1 a3 

1.3x10 


2.3x10 


^ a23 

2.7x10 


Y 


3 


345.006 


1.795 


32 


0.5 


1 1 v j 1 a3 

1.1x10 


3.9x10 


1 £.,1 a23 

3.6x10 


Y 


4 


345.506 


0.347 


22 


0.6 


-7 t„i a2 
/.2xl0 


1 A * * 1 a5 

1.4x10 


1 /z * * 1 a23 

1.6x10 


Y 


5 


345.193 


1.026 


1 1 


0.3 


3.5xl0 2 


3.3x10 s 


2.3xl0 23 


Y 


6 


345.385 


1.428 


9.8 


0.3 


3.3xl0 2 


3.3x10 s 


2.2xl0 23 


Y 


7 


345.006 


1.820 


12 


A i 

0.4 


A 1 v j 1 a2 

4.1x10 


A * * 1 a5 

2.4x10 


1 a w 1 a23 

1.9x10 


Y 


8 


345.179 


1.043 


9. 1 


0.5 


3.0xl0 2 


8.2xl0 4 


a22 

8.4x10 


Y 


9 


345.213 


1.022 


9.4 


0.4 


3.1xl0 2 


1.6xl0 5 


1.3xl0 23 


Y 


10 


344.992 


1.821 


10 


A C 

0.5 


3.4x10 


r\ - , 1 a4 

8.0x10 


*c v ✓ 1 a22 

8.6x10 


Y 


1 1 


345.115 


1.592 


6.6 


0.5 


2.2xl0 2 


7.4xl0 4 


7.0x10 


Y 


12 


345.716 


0.818 


7.2 


0.5 


2.4xl0 2 


7.5xl0 4 


7.3xl0 22 


Y 


13 


345.337 


1.018 


7.2 


A C 

0.5 


t a - 1 1 a2 

2.4x10 


1 , 1 a4 

8.3x10 


-7 w 1 a22 

7.8x10 


Y 


1 4 

14 


344.976 


1.783 


5.7 


A 1 

0.4 


1 a* * 1 a2 

1.9x10 


A -1 y. s 1 A4 

9.3x10 


-7 * * 1 a22 

7.8x10 


N 




344.980 


1.752 


3.5 


A 1 

0.3 


1 w 1 a2 
1.2X10 


A Cw1 a4 

9.6X10 


w 1 a22 
6.0X10 


IN 


1 r 


345.367 


1.434 


3.0 


A 1 

0.3 


1 w 1 a2 
1.2X10 


1 a5 
1.3X10 


c w 1 a22 
0.5X10 


Y 


1 "7 
1 / 


345.296 


1.452 


10 


A 

O.o 


1 A w 1 a2 

3.4X10 


A A vv 1 a4 

4.4X10 


c Ovyi a22 
5.OX10 


Y 


1 O 

18 


345.401 


1.419 


4. / 


A 1 

0.4 


1 /z * * 1 a2 

1.6x10 


1 1 4 ,i a5 

1.1x10 


1 wi a22 

8.1x10 


N 


19 


345.418 


1.401 


2.9 


0.3 


95 


1.9xl0 5 


l.OxlO 23 


N 


20 


345.395 


1.380 


4.5 


0.5 


1.5xl0 2 


4.0xl0 4 


4.1xl0 22 


Y 


21 


345.411 


1.406 


2.9 


0.3 


96 


1.9xl0 5 


l.OxlO 23 


Y 


22 


345.434 


1.405 


2.8 


0.3 


94 


l.lxlO 5 


7.0xl0 22 


N 


23 


345.443 


1.387 


3.5 


0.4 


1.2xl0 2 


7.9xl0 4 


6.0xl0 22 


N 


24 


345.453 


1.386 


3.2 


0.4 


l.lxlO 2 


7.4xl0 4 


5.5xl0 22 


N 


25 


345.406 


1.394 


2.9 


0.3 


95 


1.3xl0 5 


7.6xl0 22 


Y 


26 


345.389 


1.531 


2.0 


0.3 


66 


1.2xl0 5 


r c 1 A?? 

6.5x10 


Y 


27 


345.554 


1.508 


1.7 


0.2 


58 


1.4xl0 5 


7.0xl0 22 


N 


28 


345.398 


1.538 


1.7 


0.2 


55 


1.4xl0 5 


6.7xl0 22 


N 


29 


345.244 


1.021 


1.8 


0.3 


58 


9.9xl0 4 


5.5xl0 22 


N 


30 


345.301 


1.040 


1.8 


0.3 


60 


1.2xl0 5 


6.3xl0 22 


Y 


31 


345.525 


1.567 


1.8 


0.3 


59 


6.4x1 4 


4.1xl0 22 


N 


32 


345.396 


1.389 


1.5 


— 


49 


— 


~ 


Y 


33 


345.251 


1.038 


2.1 


0.3 


68 


l.lxlO 5 


6.3xl0 22 


N 


34 


344.878 


1.437 


3.y 


0.5 


1.3X1U 


A 1 1 cA 
4.1X1U 


/I 1 a22 
4.0X10 


v 
I 


35 


345.353 


1.448 


2.5 


0.5 


82 


3.1xl0 4 


2.8X10 22 


Y 


36 


345.236 


1.040 


2.1 


0.4 


70 


5.3xl0 4 


3.9xl0 22 


N 


37 


344.946 


1.796 


1.6 


0.4 


54 


3.7xl0 4 


2.8xl0 22 


N 


38 


345.442 


1.558 


1.2 


0.2 


41 


l.lxlO 5 


5.3xl0 22 


N 


39 


345.257 


1.073 


1.6 


0.3 


54 


7.4xl0 4 


4.4xl0 22 


Y 


40 


345.368 


1.036 


1.6 


0.5 


53 


1.9xl0 4 


1.8xl0 22 


Y 


41 


344.942 


1.241 


0.55 




18 






N 


42 


345.539 


1.567 


1.5 


0.3 


49 


4.5 xlO 4 


3.1xl0 22 


N 


43 


345.254 


1.052 


1.9 


0.4 


64 


5.2xl0 4 


3.7xl0 22 


N 


44 


344.904 


1.803 


1.9 


0.4 


62 


4.3xl0 4 


3.3xl0 22 


N 


45 


345.502 


0.842 


0.88 


0.3 


29 


4.4x1 4 


2.6x1 22 


N 


46 


345.559 


1.531 


0.85 


0.3 


28 


5.2xl0 4 


2.8xl0 22 


N 
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Name 


Galactic coord. 


S 1.2mm 


Diameter 


Mass 


Density 


Column density 


Infrared 




long. 


lat. 


Jy 


pc 


M 


cm -3 


cirT 2 


counterpart 


47 


345.360 


1.389 


0.80 


0.2 


26 


6.2xl0 4 


3.1xl0 22 


Y 


48 


345.441 


0.206 


1.8 


0.5 


61 


2.2xl0 4 


2.0x1 22 


Y 


49 


345.137 


1.564 


0.61 


0.2 


20 


5.7xl0 4 


2.7xl0 22 


Y 


50 


345.312 


1.046 


0.95 


0.3 


31 


5.0xl0 4 


2.9xl0 22 


Y 


51 


345.445 


1.371 


1.3 


0.5 


42 


1.5xl0 4 


1.4xl0 22 


N 


52 


345.334 


1.428 


1.5 


0.4 


49 


3.6xl0 4 


2.6xl0 22 


Y 


53 


345.563 


1.486 


1.1 


0.3 


35 


6.3xl0 4 


3.5xl0 22 


N 


54 


344.929 


1.836 


0.50 


— 


17 


- 


- 


N 


55 


345.335 


1.434 


1.1 


0.5 


36 


l.lxlO 4 


l.lxlO 22 


Y 


56 


345.434 


1.445 


1.00 


0.3 


33 


7.5xl0 4 


3.8xl0 22 


N 


57 


345.211 


1.049 


0.55 


0.3 


18 


3.7xl0 4 


2.0xl0 22 


N 


58 


345.968 


0.598 


1.5 


0.4 


49 


3.0xl0 4 


2.3xl0 22 


N 


59 


344.945 


1.229 


0.43 


- 


14 


- 


- 


N 


60 


345.006 


1.532 


1.1 


0.4 


36 


2.5xl0 4 


1.9xl0 22 


Y 


61 


345.265 


1.085 


1.2 


0.3 


39 


3.6xl0 4 


2.5xl0 22 


N 


62 


344.952 


1.195 


0.88 


0.3 


29 


6.6x1 4 


3.3xl0 22 


N 


63 


345.337 


1.037 


0.99 


0.4 


33 


1.7xl0 4 


1.4xl0 22 


Y 


64 


345.436 


1.418 


0.67 


0.3 


22 


4.0x1 4 


2.2xl0 22 


N 


65 


345.501 


0.418 


0.66 


0.2 


22 


5.4xl0 4 


2.7xl0 22 


N 


66 


345.457 


0.428 


0.81 


0.3 


27 


3.6xl0 4 


2.2xl0 22 


N 


67 


344.933 


1.251 


0.73 


0.2 


24 


7.1xl0 4 


3.3xl0 22 


N 


68 


345.064 


1.747 


0.57 


0.2 


19 


6.0x1 4 


2.7xl0 22 


N 


69 


345.501 


0.389 


1.3 


0.5 


42 


1.2xl0 4 


1.2xl0 22 


N 


70 


345.468 


1.431 


0.33 





11 


_ 


_ 


Y 


71 


345.558 


1.522 


0.40 


— 


13 


— 


— 


N 


72 


345.010 


1.764 


0.60 


0.3 


20 


1.7xl0 4 


1.2xl0 22 


N 


73 


345.590 


1.491 


0.19 


_ 


6.2 


_ 


_ 


N 


74 


344.932 


1.231 


0.23 


— 


7.5 


- 


— 


N 


75 


344.927 


1.806 


0.82 


0.3 


27 


3.0xl0 4 


1.9xl0 22 


N 


76 


344.947 


1.813 


0.48 


0.4 


16 


l.lxlO 4 


8.3xl0 21 


N 


77 


344.955 


1.170 


0.72 


0.2 


24 


6.2xl0 4 


3.0xl0 22 


N 


78 


344.936 


1.245 


0.30 


_ 


10 


_ 


_ 


N 


79 


345.384 


1.038 


0.68 


0.3 


23 


3.5xl0 4 


2.0xl0 22 


N 


80 


345.548 


1.480 


0.34 


— 


11 


— 


— 


Y 


81 


345.513 


0.407 


0.82 


0.3 


27 


4.6xl0 4 


2.6xl0 22 


N 


82 


345.563 


1.526 


0.36 


0.2 


12 


3.2xl0 4 


1.5xl0 22 


N 


83 


345.450 


1.364 


0.34 


- 


11 


- 


- 


N 


84 


345.429 


1.455 


0.88 


0.3 


29 


3.7xl0 4 


2.3xl0 22 


N 


85 


345.010 


1.617 


0.45 


— 


15 


— 


— 


N 


86 


345.001 


1.615 


1.2 


0.5 


41 


1.3xl0 4 


1.3xl0 22 


N 


87 


345.319 


1.484 


0.68 


0.4 


22 


1.8xl0 4 


1.3xl0 22 


Y 


88 


345.519 


1.639 


0.71 


0.3 


24 


4.9xl0 4 


2.5xl0 22 


N 


89 


345.133 


1.069 


1.3 


0.6 


42 


7.0xl0 3 


8.4xl0 21 


N 


90 


345.396 


1.519 


0.42 




14 


_ 


_ 


Y 


91 


345.033 


1.632 


0.76 


0.3 


25 


3.7xl0 4 


2.2xl0 22 


N 


92 


345.476 


1.568 


0.31 




10 




_ 


N 


93 


344.966 


1.181 


0.41 


— 


14 


- 


- 


N 


94 


345.240 


0.390 


0.49 




16 






N 


95 


345.524 


0.404 


0.46 


0.2 


15 


3.8xl0 4 


1.9xl0 22 


N 


96 


345.590 


0.374 


0.40 




13 






Y 


97 


345.031 


1.781 


0.56 


0.4 


19 


9.6xl0 3 


8.0xl0 21 


N 


98 


345.451 


0.435 


0.39 




13 






N 


99 


345.078 


1.786 


0.64 


0.3 


21 


1.9xl0 4 


1.3xl0 22 


N 


100 


345.854 


1.415 


0.47 




16 






N 
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Name 


Galactic coord. 


S 1.2mm 


Diameter 


Mass 


Density 


Column density 


Infrared 




long. 


lat. 


Jy 


pc 


M 


cm -3 


cirT 2 


counterpart 


101 


345.217 


0.999 


0.51 


0.3 


17 


1.8xl0 4 


1.2X10 22 


N 


102 


344.936 


1.596 


0.44 


_ 


15 


_ 


_ 


N 


103 


345.065 


1.615 


1.3 


0.5 


44 


1.7xl0 4 


1.5xl0 22 


N 


104 


345.338 


1.463 


0.73 


0.5 


24 


8.8xl0 3 


8.2xl0 21 


N 


105 


345.585 


1.484 


0.41 


0.2 


14 


3.9xl0 4 


1.8xl0 22 


N 


106 


345.854 


1.421 


0.29 




9.5 






N 


107 


345.062 


1.842 


0.53 


0.2 


18 


5.3xl0 4 


2.4x1 22 


N 


108 


345.477 


1.563 


0.22 




7.2 






N 


109 


345.454 


1.358 


0.16 




5.4 






N 


110 


345.959 


0.608 


0.33 




11 


_ 


_ 


N 


111 


344.949 


1.213 


0.27 


_ 


8.8 




_ 


N 


112 


344.942 


1.589 


0.45 


_ 


15 


_ 


_ 


N 


113 


345.539 


1.488 


0.29 


0.3 


9.5 


1.2xl0 4 


7.4xl0 21 


N 


114 


345.576 


0.258 


0.26 





8.5 


_ 





N 


115 


345.681 


0.318 


0.43 


_ 


14 


_ 


_ 


N 


116 


345.330 


1.049 


0.45 


0.2 


15 


3.9X10 4 


1.9xl0 22 


N 


117 


345.776 


1.445 


0.34 


_ 


11 


_ 


_ 


N 


118 


345.511 


1.578 


0.40 


0.3 


13 


1.6xl0 4 


9.8xl0 21 


N 


119 


344.909 


1.182 


0.65 


0.3 


22 


2.6xl0 4 


1.6xl0 22 


N 


120 


345.518 


0.835 


0.34 


_ 


11 


_ 


_ 


N 


121 


345.873 


0.811 


0.30 


_ 


10 


_ 


_ 


N 


122 


345.489 


0.430 


0.23 


_ 


7.7 


_ 


_ 


N 


123 


344.951 


1.096 


0.37 


0.2 


12 


3.7xl0 4 


1.7xl0 22 


N 


124 


344.906 


1.196 


0.38 


0.2 


13 


4.5 xlO 4 


2.0xl0 22 


N 


125 


345.092 


1.733 


0.36 




12 






N 


126 


345.261 


1.108 


0.42 


0.3 


14 


2.1xl0 4 


1.2xl0 22 


N 


127 


345.467 


0.419 


0.45 


0.2 


15 


4.8xl0 4 


2.2xl0 22 


N 


128 


344.801 


1.082 


0.27 




8.9 






N 


129 


345.482 


1.233 


0.23 




7.5 






N 


130 


345.464 


1.443 


0.21 




7.0 






Y 


131 


344.849 


1.312 


0.24 




7.9 






N 


132 


345.484 


0.856 


0.17 




5.5 






N 


133 


345.304 


0.435 


0.12 


_ 


3.9 


_ 


_ 


N 


134 


344.945 


1.600 


0.64 


0.3 


21 


2.1xl0 4 


1.4xl0 22 


N 


135 


344.950 


1.587 


0.22 




7.3 






N 


136 


345.272 


1.051 


0.36 


0.3 


12 


1.7xl0 4 


9.8xl0 21 


N 


137 


345.060 


1.727 


0.37 




12 






N 


138 


345.479 


0.400 


0.26 





8.7 


_ 




N 


139 


344.846 


1.293 


0.21 


— 


7.1 


— 


_ 


N 


140 


345.909 


0.514 


0.36 


0.2 


12 


3.4xl0 4 


1.6xl0 22 


N 


141 


345.812 


1.756 


0.26 


_ 


8.5 


_ 


_ 


N 


142 


344.957 


1.188 


0.24 


— 


7.9 


- 


— 


N 


143 


345.320 


1.417 


0.29 


0.4 


9.6 


7.6xl0 3 


5.5xl0 21 


Y 


144 


345.869 


1.820 


0.21 


_ 


6.8 


_ 


_ 


N 


145 


345.278 


1.112 


0.35 


0.2 


12 


2.7xl0 4 


1.4xl0 22 


N 


146 


345.634 


0.770 


0.33 


0.2 


11 


3.5xl0 4 


1.6xl0 22 


N 


147 


345.044 


1.618 


0.32 


0.3 


11 


2.4x1 4 


1.2xl0 22 


N 


148 


344.809 


1.075 


0.12 




3.8 






N 


149 


345.683 


1.459 


0.30 


0.3 


9.9 


2.2xl0 4 


l.lxlO 22 


N 


150 


345.566 


0.267 


0.17 




5.7 






N 


151 


344.925 


1.316 


0.33 




11 






N 


152 


345.699 


1.508 


0.18 




5.9 






N 


153 


345.286 


0.412 


0.14 




4.6 






N 


154 


345.611 


0.365 


0.27 




8.8 






N 


155 


344.916 


1.826 


0.17 




5.5 






N 
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Name Galactic coord. Si. 2mm Diameter Mass Density Column density Infrared 
long. lat. Jy pc M G cirT 3 ctrT 2 counterpart 



156 


345.472 


0.412 


0.22 




7.2 






N 


157 


345.038 


1.848 


0.28 


0.2 


9.2 


2.9xl0 4 


1.3xl0 22 


N 


158 


345.792 


1.757 


0.17 




5.7 






N 


159 


345.254 


1.012 


0.089 




3.0 






N 


160 


345.708 


1.456 


0.24 




8.0 


_ 


_ 


N 


161 


345.472 


1.231 


0.71 


0.4 


24 


1.3xl0 4 


l.lxlO 22 


N 


162 


345.497 


1.437 


0.19 


0.3 


6.4 


8.0xl0 3 


5.0xl0 21 


N 


163 


344.783 


1.105 


0.28 




9.3 






N 


164 


345.638 


0.330 


0.11 




3.7 






N 


165 


345.699 


1.853 


0.21 




7.1 






N 


166 


344.986 


1.109 


0.10 




3.4 






N 


167 


345.018 


1.762 


0.17 




5.7 






Y 


168 


344.893 


1.335 


0.45 


0.4 


15 


l.OxlO 4 


7.8xl0 21 


N 


169 


344.956 


1.348 


0.47 


0.4 


16 


6.2xl0 3 


5.6xl0 21 


N 


170 


345.948 


0.723 


0.14 




4.5 






N 


171 


344.891 


1.345 


0.20 




6.6 






N 


172 


344.924 


1.824 


0.19 




6.4 






N 


173 


344.929 


1.610 


0.15 




5.0 






N 


174 


345.247 


0.385 


0.18 




5.9 






N 


175 


345.547 


0.960 


0.13 




4.4 






N 


176 


344.892 


1.847 


0.25 




8.3 






N 


177 


345.888 


1.414 


0.26 


0.2 


8.6 


2.3xl0 4 


l.lxlO 22 


N 


178 


345.510 


0.159 


0.095 




3.2 






N 


179 


345.418 


1.469 


0.15 


0.2 


4.8 


1.7xl0 4 


7.5xl0 21 


N 


180 


344.920 


1.419 


0.17 




5.5 






N 


181 


345.079 


1.847 


0.12 




4.1 






N 


182 


345.236 


1.000 


0.20 


_ 


6.6 


_ 


_ 


N 


183 


345.375 


1.571 


0.14 





4.7 


_ 


_ 


N 


184 


345.684 


1.516 


0.14 





4.7 


_ 


_ 


N 


185 


345.532 


1.634 


0.13 





4.3 


_ 


_ 


N 


186 


345.697 


1.459 


0.19 


_ 


6.1 


_ 


_ 


N 


187 


345.501 


1.415 


0.23 


0.3 


7.5 


1.6xl0 4 


8.2xl0 21 


N 


188 


345.816 


1.765 


0.23 





7.5 


_ 


_ 


N 


189 


345.086 


1.611 


0.32 


0.4 


11 


5.2xl0 3 


4.4xl0 21 


N 


190 


345.043 


1.838 


0.21 


0.3 


7.0 


1.4xl0 4 


7.5xl0 21 


N 


191 


345.442 


1.572 


0.11 




3.8 






N 


192 


345.468 


1.225 


0.13 




4.2 






N 


193 


344.951 


1.082 


0.34 


0.3 


11 


9.6xl0 3 


6.7xl0 21 


N 


194 


345.482 


0.431 


0.27 


0.3 


8.8 


1.9xl0 4 


9.8xl0 21 


N 


195 


344.913 


1.165 


0.13 




4.4 






N 


196 


344.933 


1.820 


0.17 




5.7 






N 


197 


345.458 


1.559 


0.13 


0.2 


4.5 


l.lxlO 4 


5.3xl0 21 


N 


198 


345.792 


1.440 


0.096 




3.2 






N 


199 


345.098 


1.726 


0.11 




3.7 






N 


200 


345.293 


1.093 


0.14 




4.5 






N 


201 


344.991 


1.150 


0.11 




3.8 






N 



